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Abstract 

We study the exact one electron propagator and spectral function of a solv- 
able model of interacting electrons due to Schulz and Shastry. The solution 
previously found for the energies and wave functions is extended to give the 
spectral functions, which turn out to be computable, interesting and non triv- 
ial. They provide one of the few examples of cases where the spectral functions 
are known asymptotically as well as exactly. 

I. INTRODUCTION 

Schulz and ShastryB have introduced a new class of gauge coupled one-dimensional (ID) 
Fermi systems that are non Fermi liquid in the sense that the momentum distribution 
function has a cusp at the Fermi momentum k,F rather than a jump as in a Fermi liquid.a This 
behavior is of the sort first found by Luttinger in the context of his study of a one dimensional 
model that is popularly known as the Luttinger model.! The model introduced by Schulz and 
Shastry (SS) is in fact intimately connected to the Luttinger model, and is best viewed as a 
reinterpretation of Luttinger's original model as a gauge theory. Particles of different species 
exert a mutual gauge potential on each other, and this is sufficient to destroy the Fermi liquid. 
This model has the added property that the charge and spin correlations are unaffected by 
the interaction, owing to the "gauge" nature of interaction. The fermionic Green's functions, 
however, are non trivial, and have characteristic singularities that are popularly known as 
the Tomonaga-Luttinger liquid behavior.il! The asymptotic long distance behavior of the 
one electron correlation function is known (see below) by one of several arguments, including 
Luttinger's original one. 

Our motivation in the present work is to compute the exact one electron Green's function 
for the SS model, utilizing our knowledge of the complete spectrum of the same, and using 
techniques familiar from Anderson's treatment of the Orthogonality Catastrophe issue in 
the X-ray edge problem.§ This is of great interest since usually one does not have access 
to the exact Green's function even in ID, and one has to be content with the asymptotic 
behavior. For interpreting experiments, such as those on photoemission, one wants to know 
more than just the asymptotics, and this possibility is realized here for the particular model 
of SS. 
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We first write down the basic lattice Fermi model in ID, outline the pseudo unitary 
transformation that eliminates the gauge interactions in favor of a twisted boundary con- 
dition. Using this transformation we formulate the problem of calculating the one electron 
Green's function. 



II. THE MODEL 

Let us write the model for two component electrons hopping and interacting via the 
Hamiltonian 

L-l 

W = ~t eX P(^°" a [%,S + % + l,a]) C j,a C i + l^ + H-C. ) (!) 

3=0 <* 

where for concreteness we have simplified the original model presented in Ref. [I]. The unitary 
transformation 

U\ = exp(z a [^i,^m,i - n mA hi d )) (2) 

l>m 

transforms the Eq. ([[]) into a simple hopping Hamiltonian with twisted boundary conditions.0 
To regain a translational invariant Hamiltonian we apply a second unitary transformation 

M^Il'exp 2 ^""-^. (3) 

1=0 L 

The combined transformation U = U2U1 commutes with T, where T is the translational 
operator which shifts one site to the right (e.g. ThjT' = The effect of U on the 

fermion operators is 

UcljA ] = e-^'e^-'cj/n exp(2za<m,, ff ) IJexp 2za ( / 

z=o z=o 

= e ta ° N *e- ia ™^cJliexp(-2iaah^) JJexp ~ 2 ^(* -JW,* ? (4) 
z=o z=o 

while the density operators are invariant, Ufij a W = fij a , and the transformed Hamiltonian 
H = UUU^ reads 

n = -iEE^sWi, + H - c -) > (5) 

3=0 <T 

where n CT = N a /L is the density of a spin fermions. Thus we see that the transformed 
hopping has a "dynamically generated" gauge field. In the eigenvalue problem 

H\$) = E\~4>) (6) 

the eigenstates \(j>) are products of noninteracting one-particle states with momenta k created 
with c\ c = L~ x / 2 ^2ie lkl c[ a operator, \<f)) = H k cr c\. a \0). The momenta are quantized as 
Lkj jCr = 27rXj j(T with Zj j(T integers. The total energy and momentum of the states is 
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a j=l a j=l 



(7) 



and the one-particle energy is 

e a {k) = —2tcos(k + 2aan & ) . (8) 

Thus we must have the eigenstates of Tt 

|0) = 17+10), (9) 

with the energy and momentum given also by Eq. (|7|). In the ground state the k states 
between the Fermi momenta kp a and kp a are filled (k$ a = ±nn a — 2aan & ). In the ther- 
modynamic limit the energy E does not depend on a and is equal with the energy of the 
noninteracting a = case. For finite size systems a enters only through the 0(1/ L) correc- 
tions. 



III. SPECTRAL FUNCTIONS 



Our goal is to calculate the spectral functions, which we define as 

A a (k,uj) = E K/l4, CT |GS)| 2 5(c - Ef+' + Eqb), 
f 

B a (k,u>) = E \(f\c k jGS)\ 2 5(u; - E GS + Ef' 1 ) . 
/ 

The local {k averaged) spectral functions are defined as 

k 

1 



(10) 

(11) 



(12) 
(13) 



We concentrate on A^(k,u), since B a (k,u) is calculated analogously. 

As we mentioned in the introduction, ID interacting fermions behave as Luttinger-liquids, 
which is characterized, among others, by the power-law behavior of correlation function for 
small energies. In our case, as we will see later, the main contribution for < a < tt comes 
from 



Cl 



[uo — Bp) 2 — u 2 (k — k\ 



("AO 2 



+ Ci 



(u — Sf) 2 — u 2 (k — k 



(1)\2 



(a/Tr) 2 



uj — Bp + u(k — 



+C 2 



(a? — Bp) 2 — u 2 (k — ki ) 



W\2 



0J 



bf — u(k — ki ) 



c 2 



(uj-bf) 2 -u 2 {k-kfy 



(q/tt-1) 2 



uj — sf + u(k — fc+ ) 



(14) 



where k^> = vim a — 2ercm_ (T are the (Fermi) momenta of the singularities, c\ and c 2 are 
constants and u is the velocity of the excitations. Usually in Luttinger-liquids the velocities 



3 



of the spin and charge excitations are different and they both appear in spectral functions. 
In our case, however, due to the gauge origin of the interaction, the spin and charge velocities 
are equal to the fermi velocity v F . The spectral function has a nonanalytical, branch cut 
structure near the fermi momenta, and finite weight appears for higher multiples of the fermi 
momenta \k\ ). The local density of states near Fermi energy reads 

A(u) « Ci (uj - e F ) 2 ^ 2 + c 2 (uj - e f ) 2 ^-V 2 , (15) 

which for the noninteracting a = reproduces the Fermi-liquid step function. 

We now consider the exact evaluation of the spectral functions. As a preliminary to the 
discussion for general a, let us note the special cases of a = and a = n, where the spectral 
functions can be calculated more or less trivially. 

(i) The a = case is nothing else but the usual tight binding Hamiltonian 

n = -tJ2(cl a c J+1>a + R.c.) (16) 

of noninteracting electrons, as e icm? = 1 in Eq. (|]). For the spectral functions we recover 
the familiar 

Af(k,u) = 5(uj + 2t cos k)Q(u-e F ), (17) 
B^(k,uj) = 5(uj + 2t cos k)Q(e F -u), (18) 

i.e. a Dirac-delta peak following the cosine-like dispersion of the free fermions. 

(ii) When a = n. the model actually corresponds to the electron- hole symmetric corre- 
lated hopping modelQ with t^A — ^bb — —t and t^B = t (the hopping amplitudes t^A, ^bb 
and tAB are defined in Ref. 0): 

H = -*£(l-2n jV 0(l-2Vi,ff)4r c i+V + H - c ' ' ( 19 ) 

The Hamiltonian ( |i~9D can be diagonalized with the help of a unitary transformation 

L 

U= t (20) 

which is simpler than U = given by Eqs. (0) and (|J), and it transforms the fermionic 
operators as 

Uc\jfr = (1 - 2n^) C t. a , UCjfl = (1 - 2h^)c ha . (21) 

so the transformed fermi operators remain "local" . Furthermore, this transformation is not 
any more restricted to the ID case. The evaluation of the matrix elements is now convenient 
for operators in site representation, and the matrix element in Eq. (|T0|) becomes 

|(/|4 jT |GS)| 2 = L\(f\cljGS)\ 2 5 k>Pf . PGS , (22) 

where the c ^ creates fermion on site 0. Next, we apply the canonical transformation to 
formulate the problem using the transformed wave function [the analog of Eq. (^)], and for 
the spectral function we get 
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Mk,uj) = Lj2\(f\(l - 2u u )^\GS)\ 2 5(uj - E? +1 + i? G s)^ p;+1 _ PGS . (23) 
/ 

The wave functions are product of the spin- up and spin-down part, the evaluation is straight- 
forward and leads to 

Mk,uj) = (l-2n l ) 2 Af\k,uj) + ^ J2 E E 5^-e l {q')+e i {q)~e^k'))5 Kql . q+kl 



T2 

u g£FS| q'^FS l fe'0FS T 



(24) 



and a similar equation gives B^(k,u). In the spectral function we can identify the following 
two distinct features: (a) a Dirac-delta contribution following the cosine-like dispersion, 
which is the reminder of the noninteracting spectral function [Eq. flT7D] suppressed by a 
factor of (1 — 2n^) 2 ; (b) a broader continuum coming from the propagator dressed with a 
single loop. As we increase the filling, the weight of the Fermi jump for zero magnetization 
(n-f = ri| = n/2) decreases as (1 — n) 2 , and will disappear at half filling, leaving us with an 
A(u>) oc uj 2 density of states [the C2 ~ (1 — n) 2 in Eq. ( |i"5|) for a = ir]. To illustrate this 
behavior, we present the evolution of the local spectral functions with density in Fig. ([I]). 

(iii) General case: Like in the previous case, in evaluating the matrix elements we use 
the site representation, given by Eq. (p2[). Next, we apply the canonical transformation to 
formulate the problem using the transformed wave functions: 

(/|4 jT |GS) = (f\c^e iah ^R\GS)e- iaN i , (25) 

where R = Yii e 2ialni - l ^ L [see Eq. (|J)]. As in the transformed basis the wave functions are 
product of the spin up and down free fermion wave functions, |GS) = |GS|)|GS;) and 
\f) = l/t)|/|)> t ne matrix element factorizes, and we get 

Mk,co) = LY, |(/ T |4, t |GS t )| 2 x Kf^e^RlGS^ 2 
f 

x5(uj - Efj + £ GSiT - E fjl + E GStl )8 kt p f ^p GBA+Pf l -p GS l . 

In the equation above the cj * creates a fermion with energy e^{k') and momentum k! ^ FS-j-, 

in which case the matrix element is |(/t|cot|GSt)| 2 = 1/L. This allows us to write the 
spectral function SIS db convolution 

Af(k,u) = j £ A\{k-k',u-e^k')), (26) 



L fc'0FS T 



with 



A\(u,k) = L£|(/ ; |e^J&|GS 1 >| 5(u - E fd + E GS ^5 k , Pfd . PGSd . (27) 

h 

The interesting and nontrivial part of the calculation comes from the (/Je ia " oa ^|GS|) ma- 
trix element. In the next and crucial step, we eliminate the e iano i . This can be easily 
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accomplished after the observation that translating the operator R a similar factor appears: 
TRT^ = e^Ki-^. So 



(/jeManoa^lGS;) = e i(2m H- p ^ +PGS ^ (/| | £| GS 4 ) . 
Next, we note that e iafl °^ = (e ia + e i2a "°^)/(l + e ia ), and we end up with 



(28) 



A P ia _|_ P K P GSA- p f.l+ 2an l) - A . 

(/ile^ifllGSi) = 1 + eto (/il^iGS;) . (29) 

To evaluate (/ 4 |i2|GS 4 ), we replace |GSj_) = FT/ c^jJO) and |/ 4 ) = ]Ti%j0). Then we 

move R to the right across c£'s so that it acts to the vacuum, R\0) = |0). However, as 
Rc\^ = c*k+(2a/L) ft > the ^ momenta are shifted by 2a/ L (this is equivalent to twisting the 
boundary conditions): 



(h\R\G$i) = (oin%^ni-,ii°) = (°in^,in< +¥ , i io>- 

3=1 »=1 J=l 



(30) 



i=i 



Here we have to calculate overlap of free fermion wave functions with different phase shifts 
due to the removal of a |-spin fermion. This problem arises e.g. in the X-ray edge problem 
(Andersons's orthogonality catastrophe!), and the one-dimensional analog was discussed in 
Ref. For the reader's convenience, we repeat here the main points. The anticommutation 
relation between the operators with different phase shifts reads 



Aij 



~k l - 



e e 2 



sin a 



sin 



1 j- n 
T L 



(31) 



The overlap of the wave functions can be further calculated as |(/ 4 |-R|GS 4 )| 2 = | det A 



K/||£|GS 4 )p 



{ C fci+2a/L; C K } ' ' ' { C 1i+2q/L' C k' N } 

i C k Nl +2a/Li C *i) ••• i C k Nl +2a/L' Ck N^ 
1 



sin 2jV + q 



sin 



sin 



LL _|_ °L 
2 L 



1 



sm 



V 2 



+ 



sm ^ — H- + 1 



This determinant is actually a Cauchy determinant and can be expressed as a product, so 
we end up with 



KiOTGSj 



sm /jv i a Uj>i sm IL>i sm 



L 2N i 



n^sin 2 (^ + f 



(32) 
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For the special a = [where A+(u,k) = L5(u)5k,o] and a = 7r cases, taking the suitable 
limits, we recover the the results of Eqs. ( |17|) and respectively. In the a = 7r case the 
phase shift equals 27r/L, which is exactly the spacing between two adjacent k value, thus 
the orthogonality catastrophe is absent. 

Following the same approach, for the photoemission part we get 

B^k,u) = j £ B\{u-e^k'),k-k'), (33) 

with 

B\(u, k) = L^\(h\e- iafl ^ti\GSi}\ 2 5(u; - E GSA + E u )6 k ,p GS ^ Pfd . (34) 
/ 

The product in Eq. (^) can be evaluated numerically and spectral functions for relatively 
large system can be obtained. The numerical result is presented in Fig. ^| for some large 
size systems. Starting from a = 0, we observe that there is an overall shift in momentum 
proportional to — 2cm ^ (which we compensated for in the figure), and that apart of the main 
contribution, which follows the cosine-like dispersion, additional continuum-like features 
appear. Finally, for even larger values of a another cosine-like feature appears with a 
considerable weight. 

Alternatively, for the low energy part further analytical considerations can be applied.! 
Starting from Eq. (P2"|), the weights of the peaks can be expressed via T functions in the 
L — > oo limit, leading to the power-law behavior of the Luttinger liquid spectral function, 
and the exponents can be associated with the phase shift. We find singularities where the 
momenta of the final state are closely packed. These happen at 

k^ = unn^ — 2an^ (35) 

with v an odd integer. The most important ones for small a are with v — ±1, which 
coincides with the Fermi momenta kp+ . As we can follow in Fig. increasing a we get 

weight for the tower at k^ 3 \ which eventually becomes symmetric with k^ for a = tt, while 

the weight of the tower at /cj _1 ' ) disappears at the same time. The behavior of the primed 
spectral functions in Eq. ( p6|) has a simple behavior near k = 0: 

A\(k, u) oc [(u - e F ) 2 - u 2 k 2 ] {aM2 ~ l , (36) 

while near k = 27m j_: 

A\{k, w) oc [(w - e F f - u 2 {k - 27m;) 2 ] {a/n ' 1)2 ' 1 _ (37) 

This leads to the power law behavior of the A^(k,u) as presented in Eq. (|i"4l) . The values 
of the exponents are tabulated for some selected a in Table. |. 

The weight transfer can be quantified by observing the sum rules. While the zeroth 
momentum is constant, 



7 



B^{uj)duo = n-j- 

-oo 



A^(u)dw = 1 -n T , (38) 
the first already shows the large weight transferred to energies far from the Fermi energy: 
>B 1 (u)du = Y,(GS\cliH,c h1 ]\GS) 



UJ1 

-oo 



2t At 

— sin(7rnt) nt sin(7mjj(l — cos a), 

IX IX 



+00 



Mu)du = J2(GS\c^[H,cl 1 ]\GS) 



= — sin(7rn|) H (1 — nt) sin(7mj)(l — cos a) . (39) 

7T 7T 

We can see that the weight transfer to higher energies is the largest for a = tx and at 
half-filling. 



IV. CONCLUSIONS 

We have presented the exact one electron Green's function for a model fermi system in 
ID with a non Fermi liquid behavior for essentially any value of the interaction strength. 
A few exact analytical calculations of the spectral function for model systems, such as the 
1/r 2 exchange t-J model, with a projection to single occupancy are available in literature. i0 
The Green's function for this system obtained here does require some numerics, and is not 
totally analytical. However unlike the situation in projected models, such as the t-J model , 
it satisfies the sum rules familiar from text books for weakly interacting Fermi liquids (e.g. 
the complete electron sum rule with large u behavior of G as 1/ uj). This feature makes the 
present model particularly interesting in the context of the programme of reconstruction of 
the spectral function from its moments (e.g. see Ref. |IT|). 
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FIG. 2. The evolution of the u> and dependent spectral function as a function of a for 
n = 1/2 (left) and = 1 (right plots). The shading is proportional to A(k,uj) and B(k,uj), the 
dashed line denotes the Fermi energy. The shift of the Fermi momenta [Eq. (|35|)1 is compensated 
for by introducing q = k + an in the plot. We omitted the trivial a = case. 
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FIG. 3. The spectral function for a = 3tt/A and n = 1. Here the Fermi energy is at w = 0. 




FIG. 4. The local spectral functions B(lo) (darker) and A(u) (lighter shading) for n = 2/3. 
The a changes from (noninteracting case, top curve) to tt (bottom plot) in increments of it/ 4. 
To minimize finite size effects, the curves show the average of L = 303, 279, 255, 231, 207 and 183. 
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-6t -4t -2t e F 2t 4t 6t 

to 

FIG. 5. To illustrate the weight transfer for small a, we compare the local spectral function for 
a = 7r/4 (solid line) to the a = case (dashed). The a = 7r/4 case behaves as A(uj) ~ \uj — £f| 1//8 
near the Fermi energy. 
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TABLES 



a 





vr/4 


tt/2 


3vr/4 


7T 







1/8 


1/2 


9/8 


2 


2(a/vr- l) 2 


2 


9/8 


1/2 


1/8 






TABLE I. The exponents in the local spectral function [Eq. Jl5|)'|. 
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